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Abstract 

Fluctuations in the cosmic neutrino background are known to produce a phase shift in the acous¬ 
tic peaks of the cosmic microwave background. It is through the sensitivity to this effect that the 
recent CMB data has provided a robust detection of free-streaming neutrinos. In this paper, we 
revisit the phase shift of the CMB anisotropy spectrum as a probe of new physics. The phase shift 
is particularly interesting because its physical origin is strongly constrained by the analytic prop¬ 
erties of the Green’s function of the gravitational potential. For adiabatic fluctuations, a phase 
shift requires modes that propagate faster than the speed of fluctuations in the photon-baryon 
plasma. This possibility is realized by free-streaming relativistic particles, such as neutrinos or 
other forms of dark radiation. Alternatively, a phase shift can arise from isocurvature fluctua¬ 
tions. We present simple models to illustrate each of these effects. We then provide observational 
constraints from the Planck temperature and polarization data on additional forms of radiation. 
We also forecast the capabilities of future CMB Stage IV experiments. Whenever possible, we 
give analytic interpretations of our results. 
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1 Introduction 


Cosmology is a sensitive probe of physics beyond the Standard Model (BSM). The temperature in 
the early universe was high enough to make the production of weakly interacting and/or massive 
particles efficient. If the energy density carried by these particles was significant, then even their 
gravitational influence can be detected. This sensitivity to extremely weakly interacting particles 
is a unique advantage of cosmological probes of BSM physics. 

Until recently, the strongest cosmological constraints on BSM physics came from Big Bang 
nucleosynthesis (BBN). Measurements of the primordial abundances place stringent constraints 
both on additional forms of radiation and on injections of energy around the time of BBN (see 
e.g. [1-4] for reviews). These constraints are now complemented by the precise observations of 
the cosmic microwave background (CMB). Measurements of the temperature anisotropies have 
been used to test a wide variety of models (most recently by the Planck collaboration [5]) ranging 
from additional sources of radiation, to dark matter annihilation [6] and neutrino masses [7, 8]. In 
the future, CMB polarization experiments will play an increasingly important role in the search 
for BSM physics. Forecasts for a CMB Stage IV (CMB-S4) experiment [9] suggest an increase 
in sensitivity by one or two orders of magnitude within the next decade. With such remarkable 
improvements in sensitivity on the horizon, we should re-assess how this data could sharpen our 
understanding of the early universe, and, particularly, how it will inform our view of extensions 
of the standard models of particle physics and cosmology. 

In this paper, we will revisit the analytic treatment of the CMB anisotropies with an eye 
towards BSM applications. While numerical codes are ultimately necessary in order to make 
precise predictions for any particular model, analytic results still play a vital role. It is through 
the physical understanding of the data that we can devise new tests and motivate new models. 
For example, the use of B-modes in the search for primordial gravitational waves arose from a 
clear analytic understanding of CMB polarization [10-13]. Similarly, we wish to identify CMB 
observables that are sensitive to well-motivated forms of BSM physics, but are not strongly de¬ 
generate with other cosmological parameters. We will advocate the phase shift ^ of the acoustic 
peaks of the CMB spectrum [17] as an observable with the desired characteristics. The phys¬ 
ical conditions that lead to a phase shift are rather restrictive and determined by the analytic 
properties of the Green’s function of the gravitational potential. For adiabatic fluctuations, a 
phase shift requires fluctuations that travel faster than the sound speed of the photon-baryon 
fluid. This arises naturally for free-streaming relativistic particles, such as neutrinos (but also 
for other forms of “dark radiation”). Alternatively, a phase shift can also arise from isocurvature 
fluctuations. The phase shift therefore probes an interesting regime in the parameter space of 
BSM models (see Fig. 1). 

While particle physics experiments give strong constraints on specific scenarios, they can be 
blind to unknown or incompletely specified forms of new physics. In contrast, cosmological 
observations can provide broad constraints on phenomenological parameterizations. This has the 

^This phase shift refers to a constant shift in the locations of the high-f acoustic peaks. We emphasize that 
this is a distinct effect from the locations of the first few acoustic peaks which are sensitive to many cosmological 
parameters, as studied e.g. in [14-16]. 
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Figure 1. Particles beyond the Standard Model can be classified according to their masses M and their 
mean free paths A^fp (both normalized relative to the Hubble rate at recombination, Hrec)- Particles with 
M > Hrec contribute to the cold dark matter of the universe, while particles with M < Hrec are relativistic 
at recombination. Massive, strongly interacting particles are Boltzmann-suppressed and, therefore, do not 
contribute a cosmologically interesting density. Dark radiation separates into free-streaming and non-free- 
streaming particles. Note that axions, and other non-thermal relics, escape the simple characterization of 
this figure. 


advantage of compressing large classes of BSM physics into broad categories and is less sensitive 
to the details of the models. (This is analogous to the use of simplified models to search for new 
physics at the LHC [18], where one reduces large numbers of models to a single model which 
captures their essential features.) This approach has led to important discoveries in the past: by 
comparing observations against simple phenomenological parameterizations, the existence of dark 
matter (D^) and dark energy (Da) was established, the baryon asymmetry (r/) was identified, 
and evidence for cosmological inflation (n^) was presented. 

A useful way of parameterizing the effects of extra light species on the CMB is in terms of the 
effective number of free-streaming species, X, and non-free-streaming species, Y. The former is 
conventionally defined as the effective number of neutrinos,^ 

IVeff = a,^, (1.1) 

where Uy = ^ Here, px includes ordinary SM neutrinos (with Ny = 3.046 [20]), but 

also characterizes any free-streaming radiation density beyond the SM expectation (including 
additional energy in neutrinos). The current constraint from the Planck satellite, A"eg = 3.15 ± 
0.23 [5], represents a highly significant detection of the cosmic neutrino background. Furthermore, 
these measurements put strong limits on many extensions of the Standard Model containing 

^Stated imprecisely, colliders also constrain the number of neutrino species through precision measurements of 
the width of the Z decay. Yet, when stated more carefully, collider measurements only tell us how many fermions 
with mass below couple to the Z boson [19]. 
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additional light fields and/or thermal histories that enhance or dilute the energy in neutrinos [21- 
30].^ In addition, we will allow for a contribution from non-free-streaming radiation. We capture 
this by the following parameter:^ 

^fiuid = Oy — . (1.2) 

P-y 

We characterize the influence of A^eff and A^fluid on the photon-baryon fluid, and study their 
distinct imprints in the CMB. We also describe the types of BSM models that are being probed 
by these parameters. 

Until recently, CMB observations were not sensitive enough to distinguish between A^efr and 
A^fluid- Both types of species contribute equally to the background density of the universe and, 
therefore, affect the CMB damping tail in the same way [37]. To separate A^eff and A^fiuid requires 
measuring subtle differences in the evolution of perturbations. Free-streaming particles (like 
neutrinos) create signihcant anisotropic stress which induces a characteristic phase shift in the 
CMB anisotropies [17]. This phase shift has recently been detected for the hrst time [38]. As we 
will show, non-free-streaming particles (e.g. [30, 39-42]), in general, do not produce a phase shift 
(at least as long as the fluctuations are adiabatic and their sound speed is not larger than that 
of the photons). 

Guided by our analytic understanding, we will explore the sensitivity to these effects with the 
Planck satellite and with a future CMB-S4 experiment, focusing on the ability to distinguish the 
parameters A'eff and A^fluid- Our analysis of the Planck temperature and polarization data leads 
to the following constraints:® 

A^eff = 2.801°;24 (la), A^fl.id < 0.67 (2a). (1.3) 

We see that the current data is already sensitive to the free-streaming nature of the fluctuations. 
We will explain the important role played by the polarization data in breaking the degeneracy 
between A'eff and A^fluid, as well as that with the helium fraction Yp. We will also show that a 
CMB-S4 experiment would improve these constraints by up to an order of magnitude under a 
number of experimental conhgurations. We will highlight how present and future constraints are 
driven by measurements of the phase shift. We will also explore how these measurements may 
be optimized, including through the use of delensing to sharpen the acoustic peaks. 

The outline of the paper is as follows. In Section 2, we derive analytically the effects of new 
relativistic particles on the perturbations of the photon density. We identify the precise physical 

^For thermally-decoupled massless particles, the equivalent change to N^b will be a function of the decoupling 
temperature. A single thermally-decoupled species implies SN^b = N^b — 3.046 > 0.02 [25]. Interestingly, this level 
is within reach of future CMB-S4 experiments [9]. 

^Another attempt to parametrize free-streaming radiation is in terms of a viscosity parameter Cvis [31]. This 
parameter has recently been detected by Planck, Cvia = 0.331 ± 0.037 [5] (see also [32, 33]). However, as discussed 
in [34], cJia = I is not equivalent to free-streaming radiation and differs from ACDM by = 20. Our pa¬ 
rameterization has the advantage that it reproduces ACDM when N^b = = 3.046 and Afiuid = 0. A similar 

parameterization has appeared in [30, 35, 36[, and was analyzed with WMAP data in [35, 36]. However, it has 
only recently become possible to distinguish these parameters with high significance. 

®The constraints assume that the helium fraction Yp is fixed by consistency with BBN. Results that marginalize 
over Yp are presented in §3.2. 
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conditions that produce a phase shift in the CMB anisotropy spectrum. We illustrate these 
effects through an exactly solvable toy model. We also compute the phase shift for a simple 
model with isocurvature fluctuations and for free-streaming relativistic particles. In Section 3, 
we confirm some of these pen-and-paper results through a numerical analysis. We present new 
constraints on dark radiation from the Planck 2015 data [43] and forecast the capabilities of 
future CMB-S4 experiments [9]. Section 4 contains our conclusions and a description of plans 
for future work. In two appendices, we comment on the inclusion of matter (Appendix A) and 
polarization (Appendix B) in our analytic treatment. 

Notation and Conventions 

We will work in conformal Newtonian gauge 

ds^ = a^(r) [(—1 — 2‘h)dr^ + (1 “ 2'I')(5jjdxMx'^] , (1.4) 

where r is conformal time. We use tq for the present time, Tree for the time of recombination, Teq 
for matter-radiation equality and Tin for the time at which we set the initial conditions. We will 
use a subscript ‘a’ to denote quantities evaluated at the time Ta- The conformal Hubble parameter 
18 % = a/a, where overdots stand for derivatives with respect to t. It will be convenient to define 
the sum and the difference of the two metric potentials, 


$± = $±4^. (1.5) 

We will use p and P for density and pressure, and a for the scalar potential of the anisotropic 
stress. Individual components (like photons, matter, neutrinos, etc.) will be denoted by a sub¬ 
script a = 7 , m, I/, • • •, where X and Y will stand for free-streaming and non-free-streaming 
species, respectively. We use overbars for homogeneous background quantities, pa, and write 
perturbations as 6 pa = pa — pa- Following [17], we define the overdensity in the particle number 
with respect to both the proper volume, 5a = 5pa/{pa + Pa), and the coordinate volume, da- The 
relation between the two definitions is da = 5a — 34'. The equation of state and the speed of 
sound are Wa = Pa/pa and = 5Pa/5pa, respectively. 
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2 Physical Origin of the Phase Shift 

The structure of the acoustic peaks in the CMB is largely determined by the propagation of 
fluctuations in the photon-baryon plasma. The physics is that of a harmonic oscillator with a 
time-dependent gravitational forcing, 

, ( 2 . 1 ) 

where ~ A non-trivial evolution of is sourced either by anisotropic stress a or by pressure 
perturbations 6 P (see Fig. 2). Under certain conditions, which we will identify, this induces a 
contribution to d^y which is out of phase with its freely oscillating part. In this section, we will give 
an analytic description of these effects, building on the pioneering work of Bashinsky & Seljak [17] . 


pressure perturbations (§2.3, 2.4.1) 

5P 
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anisotropic stress (§2.4.2) 


Figure 2. Illustration of the coupled perturbations in the primordial plasma. 


2.1 Preliminaries 


We begin by collecting a few standard results from cosmological perturbation theory (see e.g. [17, 
44, 45] for further details). This mainly serves to fix our notation and to introduce the main 
equations used in this paper. 


The CMB couples gravitationally to perturbations in the matter fluctuations. We define the 
stress-energy tensor for each species a as 

T%,a = -{pa + Spa), T\,a = {Pa + Pa)Vi,a , Tj,a = {Pa + 6Pa)5) + {pa + Pa)^^j,a ■ ( 2 . 2 ) 


The scalar part of the velocity can be written as Vi^a = —PiUa, where Ua is the velocity potential. 
Similarly, the anisotropic stress tensor can be expressed as = ^{PiPj ~ 
where the factor of | was introduced for future convenience. Conservation of the stress-energy 
tensor for each decoupled species implies 

L + 3T , (2.3) 

Ua = cl6a - XaUa + V^CTa "b , (2.4) 
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where Xa = ^(1 — 3c^) is the Hubble drag rate. Equations (2.3) and (2.4) can be combined into 


da + Xada “ = V'^CJa + V^($ + 3c^4') . 


(2.5) 


To solve this equation requires additional equations for a a and 

Boltzmann Equation 

An evolution equation for Ua can be derived from the Boltzmann equation for the distribution 
functions fa{T,x,q,h) of each particle species a with comoving momenta q = qh. We separate 
fa into a background component fa and a perturbation 5fa = fa — fa- For massless particles, it 
will be convenient to integrate out the momentum dependence of the distribution function and 
define 



( 2 . 6 ) 


The linearized, collisionless Boltzmann equation is then given by 


Da + n'ViDa = -3h*Vi4>+ . 


(2.7) 


We note that Da only depends on <h+, but not on 4>_. It will be useful to expand the momentum- 
integrated distribution function Da into multipole moments. 


OO 



( 2 . 8 ) 


where the Legendre polynomials Pi{lf) are functions of ^ = h ■ k. The monopole moment 
Dafi determines the overdensity da, while the quadrupole moment Da ,2 is associated with the 
anisotropic stress Ua. To see this, one writes the perturbed stress-energy tensor in terms of the 
perturbed distribution function. 



(2.9) 


Comparing this expression to (2.2), we find 



( 2 . 10 ) 


The quadrupole moment of (2.7) then provides the missing evolution equation for the anisotropic 
stress. 


Einstein Equations 

The metric potentials <I> and T are determined by the Einstein equations 


T + 'H{24> + 6) + {2n + + Jv2($ - ^) = ATTGa^SP , 

o 

- 3n{A> + nA>) = AttGo^ 5p , 


( 2 . 11 ) 


( 2 . 12 ) 
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where 5p = Yha ^Po- = Sa are the total density and pressure perturbations, respec¬ 

tively. In terms of the fields <I>±, the evolution equation (2.11) becomes 

$+ + m^+ + {2% + n^)^+ = SttGo^ 5P + 5[^>_], (2.13) 

where 5[<I>_] = +1-1^- — {2'H + 'H^ + |V^)<I>_. The held $_ is related to the total anisotropic 

stress a via the constraint equation 

<!>_ = —127rGa^ (p-|-.P)(T , (2.14) 


where (p + .P)cj = YlaiPaPPa)o'a- In the standard model, both photons and neutrinos contribute 
to 6 P, but only neutrinos provide a source for a (and hence $-). BSM particles may lead to 
additional pressure and/or anisotropic stress. 

During radiation domination, the evolution equation (2.13) can be written as 




4 

y 


SttGo^ 

{c-^kY 



c^) 5pa + 5[$_], 


(2.15) 


where primes denote derivatives with respect to y = c^kr and 5[<h_] = 4>'/. + (2/y)4>'_ + 34>_. 
The Green’s function for (2.15) is 


y" L 


y 

G<i>+ {v, y') = Q{y-y')^ {y' - y) cos(y' - y) - (l + yy') sin(y' - y) 


(2.16) 


where 0 is the Heaviside function. Given a model for the sources in (2.15), we use the Green’s 
function to determine the solution for <1>+. The time evolution of the source terms may require 
solving the Boltzmann equation (2.7). 


2.2 Conditions for a Phase Shift 


We are now in the position to analyze the evolution of perturbations in the photon-baryon plasma. 
For simplicity, we will ignore the small effect due to the baryons,® so that the speed of fluctuations 
in the photon density is The Hubble drag rate in (2.5) therefore vanishes, X 7 ~ 0, and 

we get the evolution equation (2.1). The solution for can then be written as 


cf-y(y) = in cosy - 


dy'4>+(y')sin(y-y') 


(2.17) 


where the hrst term is the homogeneous solution with constant superhorizon initial condi¬ 
tion = d^{y\n ^ !)• The second term is the inhomogeneous correction induced by the 
evolution of metric fluctuations. Since sin(y —y') = siny cosy' —cosy sin y', we can write (2.17) as 


d^{y) = d^^in + c A{y) cosy-c ^B(y) siny 


(2.18) 


°We ignore the contributions of baryons and dark matter to the energy density, but we are implicitly including 
the baryons when we assume that the photons are not free-streaming particles. 
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where 


My) = 

rv 

/ dy'^>+(y') siny, 

Jo 

(2.19) 

Biy)^ 

ry 

/ dy'^>+(y') cosy'. 

Jo 

(2.20) 


We wish to evaluate (2.18) at recombination, y —)• y^ec- For the high-£ modes of the CMB, it 
is a good approximation to take the limit y —)• oo and assume that the background is radiation 
dominated (see Appendix A for further discussion). If the integral in (2.20) converges in this 
limit, then a non-zero value oi B = limj^_^oo B{y) will produce a constant phase shift 9 relative 
to the homogeneous solution, 

sin(9 = ^= - . (2.21) 

(yA-\- +3"^ 


This phase shift will be reflected in a shift of the acoustic peaks of the CMB anisotropy spectrum. 
In the following, we will identify the precise physical conditions for which such a phase shift is 
generated. 

It will be convenient to combine B and A into a complex field 


B + iA 



dy ^+{y) 





( 2 . 22 ) 


where $^^(y) is an even function of y, while ^>^^(y) is an odd function. It is easy to see that 


the even part of <h+ determines B and the odd part determines A: 


B = - 
2 


1C 

2 7 -. 






f 


dy $^^(?/) 


(2.23) 


(S') ' is^ 

We will get B = 0 as long as <I>Y {y) is an analytic function and (y) vanishes faster than 

y~^ for \y\ —?> oo.' This suggests two ways of generating a non-zero B and, hence, a phase shift 
in the solution for the photon density: 


(S) 

i. rapid growth of (dziy) —?■ mode traveling faster than 

ii. non-analytic behavior of <kY(y) —non-adiabatic fluctuations. 


The mathematical requirements listed on the left are mapped directly into physical conditions, 
shown on the right. 

• The first condition is easy to understand physically: in (2.22), the Green’s function of d^, 
i.e. sin(y — y'), leads to exponential suppression for y —)> ioo. To have a growing solution at 
y = foo, we therefore need a term in of the form _ ^-i(cslc^)y 8 

^Since the equations are symmetric in i/ —> —y, the odd part analytic around j/ = 0. This is why 

we always find contributions to A. 

®Note that Cs is just a parameter of the wave-like solution and is not necessarily the sound speed of a fluid. 
Indeed, in the case of free-streaming radiation, it corresponds to the propagation speed of the individual particles. 
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• The second possibility, non-analyticity, is easy to understand mathematically, but the phys¬ 
ical requirements are less transparent. First of all, the equations of motion for any mode 
should be analytic around any finite value of kr in the radiation-dominated era since there 
is no preferred time. Hence, the only moment at which non-analytic behavior is possible is 
around kr = 0, i.e. where the initial conditions are defined. Let us first show that adiabatic 
initial conditions are analytic at kr = 0. By definition, for adiabatic initial conditions, any 
long-wavelength mode is locally generated by a diffeomorphism [46]. In the limit fer —)• 0, 
we then have = 4>_|_^in-|-0(fe^r^). This expansion is necessarily analytic in (by locality 
and rotational invariance), but also in because the scaling k ^ Xk and r —)■ A“^r 

can be absorbed into the overall normalization of the scale factor a which has no physical 
effect.® Hence, must be analytic around y = c^kr = 0, as long as the modes are 

adiabatic. Conversely, any violation of analyticity requires a source of non-adiabaticity. 

In the following sections, we will illustrate the different physical origins of the CMB phase shift 
through a number of simple examples. 

2.3 Intuition from a Toy Model 

To gain more intuition for the system of equations discussed in the previous section, let us solve 
them exactly in a simple toy model. In particular, we will study an example in which the 
metric fluctuations propagate with a different speed than the photons, Cs We wish to 

understand under which conditions this mismatch leads to a phase shift in the photon oscillations. 

In the absence of anisotropic stress, the Einstein equation for <I>+, eq. (2.13), is 

<1+ + + {2n + n‘^)^+ = sp . ( 2 . 24 ) 

The Friedmann equations describing the evolution of the homogeneous background imply 

2n + n^ = -8TrGa^P = -3n^w, (2.25) 

where we have defined the equation of state w = P/p. (Note that we are not assuming that 
P = wp.) We write the pressure perturbation as 

5P = cl5p + 5Pen, (2.26) 

where Cs is the speed controlling the propagation of the total density perturbation 5p and <5Pen 
denotes the non-adiabatic entropy perturbation. For adiabatic fluctuations, one has dPen = 0 
and = w — [3?^(1 -|- rc)]“^ w. We eliminate the density perturbation 6 p using the relativistic 
generalization of the Poisson equation (2.12). Eqnation (2.24) can then be written as 

+ 3-^(1 + 0^)6+ - 3n'^{w - c2)$+ + cik^<^+ = SvrGa^ hPen ■ (2.27) 

So far, this is fairly general and has only assumed vanishing anisotropic stress. In particular, 
at this point w and Cg are still general, possibly time-dependent parameters. To be able to derive 

®In a universe with a preferred time, this rescaling would also require a shift in this preferred time to keep the 
density fluctuations fixed. For adiabatic modes, the curvature perturbation ( is conserved outside the horizon even 
in the presence of such a preferred time and so this is unlikely to have an impact on gauge-invariant observables. 
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an analytic solution for the evolution of <i>+(r), we will now make a few simplifying assumptions. 
First, we assume that the equation of state w is nearly constant, so that we can integrate (2.25) 
to get 


n 


2 1 
1 + 3w T 


(2.28) 


Second, we take Ri const, and 6 Pen ~ 0. This allows us to solve (2.27) analytically. For 
arbitrary w and Cg, these assumptions are not guaranteed to be easily realizable in a physical 
model. Our analysis only serves as a simple illustration of some of the effects that give rise 
to phase shifts in the CMB. More concrete examples of these effects will be discussed in §2.4.1 
and §2.4.2. 


It is convenient to define z = Cskr and write (2.27) as 


dz^ 


4 >+ + 


1 - 2 a 
z 


d^ 


4 >+ + 



$+ = 0 , 


(2.29) 


with 


_1 3(1+ c^) _ 12{w - eg) 

2 l + 3w ' ^ (l + 3'u;)2 ■ 


(2.30) 


In the physically interesting parameter regime, 0 < (Cg,w) < 1, we have — ^ < a < — |, with 
equality for {c^,w) = (1,0) and (0,1), respectively. The general solution of (2.29) is 


^+{z) = {ciJ^,{z)+ C2 Y^{z)) , At = + /3 , (2.31) 


where Jk(z) and Yk(z) are Bessel functions of the first and second kind, respectively. Note that At is 
strictly positive and real-valued for the physically relevant parameter regime: ^V3 < At < 
with the minimum at {c^,w) = (^, 1 ) and the maximum at {c?g,w) = ( 1 , 0 ). 

To impose initial conditions, we consider the super horizon limit, z <C 1, 


2 "^(Cl +C2COt(7rAt)) 2'=C2F(At) ^ 

F(l + At) vr ^ + 


(2.32) 


Since a -|- At > a — At, the “growing mode” solution corresponds to C 2 = 0. Hence, we have 


$+(z) = ciz“J«(z). (2.33) 

The overall normalization in (2.33) will depend on the nature of the initial conditions (adiabatic 
or entropic). For = w, the superhorizon limit of is a constant, which we match to the 
super horizon value of the primordial curvature perturbation This leads to the normalization 
Cl = 2\f^ Q. We will maintain this normalization even for c^ 7 ^ rc, although, in principle, the 
normalization of non-adiabatic modes is model-dependent. The y —)• 00 limit of (2.20) then 
becomes ^ 

B = 2V^cj dy J^^^y^cosy, (2.34) 

and similarly for A. Together with d^^in = —3^, this allows us to compute the phase shift 
via ( 2 . 21 ). A graphical illustration of the dependence of the phase shift 9 on the parameters 
cl and w is given in Fig. 3. Let us emphasize again that we do not imagine that all of the 
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Figure 3. Phase shift 6 for varying speed of sound (cs) and equation of state {w). The dashed line 
denotes 9 = 0. Below this line, the phase shift is negative, while above it is positive. 


combinations of and w that we show in the figure can be realized in a physically realistic 
model. In the following, we take slices through the parameter space to show that the most 
important features of the figure can be understood analytically. 

Consider first the special case = w, which corresponds to adiabatic fluctuations. The 
parameters in (2.30) and (2.31) then reduce to 

5 + 3c^ 

“^~2(rT3^’ ^ = K=\a\. (2.35) 

As shown in the left panel of Fig. 4, the phase shift then vanishes for Cg < c-y = ^ and is 
positive for Cg > c^y. This is consistent with our abstract reasoning of the previous section. At 
large z = Cgkr, the solution (2.33) behaves as cos(z) oc cos{cg/c-y y), with y = c-ykr. Since 
the contour at infinity in (2.23) will not vanish when Cg > Cj, we cannot conclude that 0 = 0 
(cf. condition i. in §2.2). To find 0^0 was still not guaranteed, but there was no reason to 
expect otherwise. In contrast, 0 vanishes for Cg < for exactly the reasons discussed before. In 
particular, the solution (2.33) now takes the form with a < 0. Near z = 0, the solution 

is analytic (cf. condition ii. in §2.2) since the non-analytic behavior of z“ cancels that of the 
Bessel function. Of course, this is precisely what we expected for adiabatic modes. 

Taking ^ w corresponds to non-adiabatic fluctuations, i.e. fluctuations which contain an 
isocurvature component. In this case, we expect a phase shift to arise for any values of and w. 
To Illustrate this, let us take = | and only allow to vary. We then have 

3^2 _ 

a = -l-^, P = l-3cl, K = -y/8 + 9cj. (2.36) 
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Figure 4. Phase shift 9 for varying Cg = w (left) and for varying and fixed w = ^ {right). 


The corresponding phase shift is shown in the right panel of Fig. 4. We see that the phase shift 
now only vanishes at the special point = tc = |, where the fluctuations are adiabatic. This is 
also where the parameter /3 changes sign, which is the origin of the change in the sign of the phase 
shift, 0^0 for Cg ^ tc. This time the phase shift is associated with the non-analytic behavior of 
<h+(z) near the origin. To see this explicitly, consider the z —>■ 0 limit of (2.33): 

= 2.^ra + „) 

For 0 < Cg < 1, we have a + k <2 and, hence, 4 )+( 2 ) is non-analytic at z = 0 (cf. condition ii. 
in §2.2), except for the special case of the adiabatic limit where a + k = t). This corresponds to 
the non-trivial super horizon evolution of 4>+ in the presence of isocurvature modes. 

2.4 Simple Examples 

2.4.1 Isocurvature Perturbations 

The toy model of the previous section suggests that isocurvature perturbations produce a phase 
shift. To study this further, it is useful to consider a slightly more realistic model. To simplify the 
calculations as much as possible, our curvaton-like model will include three species: photons ( 7 ), 
a dark radiation fluid {Y) and a matter component (m) that decays into the dark radiation. The 
matter will carry the isocurvature fluctuations. We include the dark radiation because we are 
only interested in the gravitational effects on the photons, i.e. we want to avoid the direct heating 
of the photons due to the decaying matter. The model will allow us to explore small deviations 
from the adiabatic limit = w. 

The coupled equations for the background densities of m and Y are 

X-^ia^Pm) = (2.38) 
dr 

-^^{a’^py) =+Tapm, (2.39) 

a* dr 
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where F is a constant decay rate. To simplify calculations, we will work perturbatively in 
Cm = Pm/p- At zeroth order in e^, the universe is radiation dominated, and hence a = tIt^. 
Integrating (2.38), we get 




(2.40) 


where we set the initial value Pm,\n = Pm{T\n)- Substituting (2.40) into (2.39), we would get the 
solution for py{cl), however, this will not be needed for our purposes. 


We now wish to determine how the decaying matter affects the evolution of the metric pertur¬ 
bation In the absence of anisotropic stress, this is given by (2.24). The pressure perturbations 
only receive contributions from 7 and Y, so we have 


5P = c^^{6p^ + 6pY) = -Sp--6pm, 


(2.41) 


where we have used Cy = = |. Using the Poisson equation (2.12), we can write (2.24) as 


$+ + - - V^4>h 

O 


(Sw - - nhmSm ■ 


(2.42) 


We wish to solve this at linear order in e^- 

We will shortcut the computation by isolating the isocurvature contribution. Suppose we write 
and similarly for 6 m- Equation (2.42) then implies 




+ Z^iso/ ^ 

y 


riso I \ 

o ^rn ^\^m) ? 


(2.43) 


where we have used that $ 1 ^° ~ 0 {em), so that all terms multiplying can be evaluated at 
zeroth order in e^. Since the right-hand side of (2.43) is proportional to e^, we only need the 
homogeneous solution for 5)^°, which is 


C°(y) = Cl + C2lny , 


(2.44) 


where C 17 are constants that may depend on k. We solve (2.43) using the Green’s function (2.16). 
Substituting (2.40) and (2.44), we get 




C-ykTm 


'Vin 


f 

[ 1 iy'? 1 

( - exp 

2 (c-yfeT^iec)^ 


Cl C 2 ln(?/') 


(2.45) 




where im,in = em,ine 2 ^’^“Adec) have introduced the “decay time scale” = rin/T. 

Let us comment on a few features of this solution. First of all, we notice that the integral 
is highly suppressed when krdec ^ 1- The reason is easy to understand: the integral would 
have been dominated by contributions around the time of horizon crossing, y ~ 0 ( 1 ), but, for 
kTdec ^ 1) this is long after pm has decayed. Second, we see that the solution has an overall factor 
of {c-ykTm)~^■ This reflects the growth of from the initial time. Tin, to the time of horizon 
crossing, {c-yk)~^. 
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It is convenient to define Tgq as the time at which pm and + py would be equal if there was 
no decay. This is given by Teq ~ Equation (2.45) then becomes 

= A-m ■ (2-46) 

C'y/CTeq 

We compute the phase shift by substituting (2.46) into (2.19) and (2.20), and taking the limits 
yin 0 and y oo, 

1 1 '°° 

— / dy'X(y')siny', (2.47) 

C-y/CTeq Jq 

= / dy'X(y')cosy'. (2.48) 

C'y/CTeq Jq 


In Fig. 5 we display the numerical result for as a function of ydec = c^A:rdec- Iii the limit 
Vdec ^ 1? we can simplify the calculation by dropping the exponential in (2.45). We then get 


^iso ?/dec^l ^ 
^iso ^dec^l ^ 


TT 1 

4^^ C^/cTeq ’ 

^(ci -C27 e) 


1 


Cq/A^Teq 


(2.49) 

(2.50) 


where 7 e ~ 0.5772 is the Euler-Mascheroni constant. We see from Fig. 5 that the analytic 
result (2.50) becomes a good approximation for ydec A 5- 



Figure 5. Numerical value of as a function of ydec- The blue and red solid lines show the 

effect from ci and C 2 , respectively. The dashed lines are the asymptotic values calculated in (2.50). 


To summarize, we have demonstrated in a simple model that isocurvature perturbations give 
rise to a phase shift, as we expected from condition ii. of §2.2. As suggested by Fig. 5 this phase 
shift has a nontrivial scale dependence which probably allows it to be distinguished from other 
sources for a phase shift. This scale dependence is likely to be a general feature of isocurvature 
models. 


14 












2.4.2 Free-Streaming Particles 


Above we have seen that a phase shift is also generated if fluctuations in the gravitational potential 
travel faster than the speed of sound in the photon-baryon fluid. A simple way to realized this 
is through free-streaming relativistic particles, such as neutrinos [17]. In this section, we give 
a new derivation of this result. In Section 3, we will show that the CMB data is now accurate 
enough to detect this effect (see also [38]). 

Since most of the modes relevant to current and future CMB observations entered the horizon 
during the era of radiation domination, our discussion in this section will ignore both the matter 
and baryon content of the universe. In Appendix A, we show that this is a good approximation 
for high-^ modes and also discuss some of the implications of a finite matter density. 

We start with the Einstein equation during radiation domination, eq. (2.15), which in the 
absence of additional pressure perturbations 6 Pa = c^5p, with Ca 7 ^ c^, takes the form 

-h -4>', -h 4>+ = 5[4>_] 

y 

= +3$_ . (2.51) 

y 

Using the Green’s function (2.16), the solution for <1>_|_ can be written as 

= r <^yS[^-{y')]G^Ay^y')- ( 2 - 52 ) 

y Jy-m 

Following [17], we will write this as an expansion in powers of the fractional energy density 
contributed by the species of free-streaming particles, X, as measured by the dimensionless ratio 

« - ■ (2.53) 

P') + Px o-u + Agff 

For the Standard Model neutrinos, we have ~ 0.41. We determine the superhorizon initial 
condition of the homogeneous solution, <h+,in, by matching to the constant superhorizon solution 
for adiabatic initial conditions [17], 


+ = 


20 + 4ex 
15 + 4ex ^ 




(2.54) 


where C is the conserved curvature perturbation. 

To find the inhomogeneous part of the solution (2.52), we require <I>_(y). This is determined 
by the anisotropic stress ax induced by the free-streaming particles, 

(y) = - crx{y) = -\ ^Dx, 2 {y) , (2.55) 

y 3 

While it should be physically clear that free-streaming radiation travels at the speed of light, this property 
is sometimes not very transparent in the equations for the density perturbations of this radiation. Instead, what 
is usually more apparent is that free-streaming particles can induce a significant anisotropic stress, which then 
provides a source for 4+ and, hence, the evolution of d~, through (2.1). The origin of the phase shift is therefore 
often identified with the presence of anisotropic stress. However, in principle, one could imagine situations with 
significant anisotropic stress, but no supersonic propagation modes (e.g. non-relativistic, free-streaming particles). 
In those cases, we would not expect a phase shift. Hence, it is the propagation speed, not the anisotropic stress 
itself, that makes the phase shift possible. 
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where Dx ,2 is the quadrupole moment of the momentum-integrated distribution function defined 
in (2.6). The first equality follows from the Einstein constraint equation (2.14) with a = exO’X; 
whereas the second equality employs (2.10). The evolution of Dx ,2 is determined from the 
linearized, collisionless Boltzmann equation (2.7), 

Dx + ikfiDx = —Mkfj. . (2.56) 


Defining Dx,m 


Dx{Tin) at some time Tin, the solution to (2.56) is 


Dx{t) 


— g-*fc/i(T-'rin) 


Dx,in — 



dr' e-*^^(^-^')$+(r'). 


(2.57) 


We wish to extract the quadruple moment Dx ,2 of the solution. Since Dx,r(Tn) oc Tj^, we will 
only keep the monopole term Dx,o(Tn) in the homogeneous part of the solution. This is possible 
because we can set the initial conditions at a sufficiently early time Tin A:”^, so that the modes 
with £ > 0 will be subdominant. In fact, we will take fcrin —)• 0 from now on. Assuming adiabatic 
initial conditions, i.e. Dx,o(Tn) = dx,m = —3^, we get 


Dx, 2 {y) = -Kj 2 [c^^y\+ ^ dy'4>+(?/') |^ji[c^ny - y )] - - 2 / 0 ] I > (2-58) 

where the Bessel functions arise from the Rayleigh expansion of the exponentials. Substituting 
this into (2.55) directly links the two gravitational potentials <1>+ and <1>_. The most important 
feature of the solution (2.58) is that it contains modes that travel at the speed of light. Specifically, 
recall that c~^y = kr and, therefore, the Bessel functions describe oscillatory solutions with a 
speed of propagation of c* = 1. As we have emphasized before, this is the property of the 
free-streaming radiation that makes a phase shift possible. 


The above is a closed set of equations which we can solve perturbatively in ex' 

= = (2.59) 

n n 

(n) (n) 

where the superscripts on <k)|_ ^ and d\ ' count the order in ex- Here, we present the solution up 
to first order; 


At zeroth order in ex-, we have = 0 and, hence, 

solution, 

^(0)(^)^4^siny-j^ 


(0) 


y-' 


Inserting this into (2.19) and (2.20), we find 


A<")fe) = 2C-2C?4!<h J£ 


i?W(y) = 2C 


2 /" 

y — cos y sin y y^oo 


y" 


The result for the photon density perturbations then is 

4 °^(y) SCcosy. 


is given by the homogeneous 



(2.60) 

2C, 

(2.61) 

0. 

(2.62) 


(2.63) 
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We conclude that in the absence of anisotropic stress, the correction due to $+ is in phase 
with the homogeneous solution and B vanishes as expected. 


• At first order in ex, we only need the zeroth-order solution of the anisotropic stress, 
since the source in (2.55) already comes with an overall factor of ex- Hence, eqs. (2.55) 
and (2.58) can be written as 


c^H)(y) 


( 2 - 64 ) 

~ ^ ^ -y')]- 5^3 [c-\y-y')]^ , 


where <i>® is given by (2.60). Substituting (2.64) into (2.52), we obtain 


A"!!/) = 


sm y — y cos y 


+ dy'5[4>L'^(y')]G'$^(y,y')- 


r JO 

Inserting this into (2.19) and (2.20), we finally get expressions for AH) and 
have to be evaluated numerically, and we find 


(2.65) 

These 


« -0.268 C ex , B^^^ ^ 0.600 C ex ■ 


( 2 . 66 ) 


The non-zero value of HH) corresponds to the expected phase shift. 


Using (2.21) with d-y^in = —3C, we get 

0«O.19l7rex + O(ei), (2.67) 

which is consistent with the result of Bashinsky & Seljak [17]. 

The phase shift is a clean signature of free-streaming particles and will naturally play an 
important role in the observational constraints discussed in Section 3. To put these constraints 
into context, let us use the analytic result of this section to relate changes in iVeff to shifts 
5i in the peaks of the CMB spectra. As we show in Appendix B, the E-mode spectrum will 
exhibit precisely the same phase shift as the temperature spectrum and, therefore, our analytic 
estimates are applicable in either case. In the flat-sky approximation, a shift in angle 9 is related 
to a multipole shift by 6 i ~ (d/vr) A-^peak, where A^peak ~ 330 [43] is the distance between 
peaks in the temperature anisotropy spectrum for modes entering the horizon during radiation 
domination. Using (2.67), with Agg = = 3.046, we find that the shift of the peaks arising from 

ordinary neutrinos is 5ii, ~ 26 (compared to a neutrinoless universe). Similarly, small variations 
in Agff around the standard value will lead to a multipole shift of order 

~ 5.0 X 5Nes , (2.68) 

where we have expanded to linear order in dAgg = Agg — 3.046. While this result is likely subject 
to a 20 to 30 percent error, it is reliable enough to see that a sensitivity of a{Nes) ~ 0.1 will 
constrain a phase shift of order <5£ < 1. Current constraints on Agg imply 5i ^ 0(1). As we 
will see in Section 3, future CMB experiments are expected to constrain, or measure, shifts of 
order 5i ~ 0(0.1). This is consistent with the rough expectation from measuring 0(10) peaks 
and troughs in the E-mode power spectrum. 
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3 BSM with the CMB: Present and Future 


The CMB has the potential to distinguish between many distinct sources of BSM physics; new 
free-streaming or non-free-streaming particles, isocurvature perturbations and/or non-standard 
thermal histories. However, ultimately the observability of the new physics depends both on the 
size of the effect and whether it is degenerate with other cosmological parameters. In this section, 
we present new constraints on the density of free-streaming and non-free-streaming radiation 
from the Planck satellite [43] and then discuss the capabilities of a proposed CMB-S4 mission [9]. 
Whenever possible, we will give some approximate analytic understanding of the qualitative origin 
of our results. For precise quantitative results, we will perform a full likelihood analysis. 


3.1 A Discussion of Degeneracies 

It is useful to anticipate the possible degeneracies between the effects of extra relativistic species 
and changes in the cosmological parameters. Here, we show that the degeneracy with the helium 
fraction Yp can be understood analytically. 

At the level of the homogeneous background, the largest effect of relativistic particles is a 
change of the expansion rate during radiation domination, 

+ , (3.1) 


where free-streaming particles (Aeff) and non-free-streaming particles (Afluid) contribute equally. 
The change in the Hubble rate manifests itself in a modification of the damping tail of the 
CMB. Specifically, for modes with wavelengths smaller than the mean free path of the photons, 
A < Amfp, the temperature fluctuations are suppressed by exjp[—{k/kd)‘^], where kd is the wave 
number associated with the mean squared diffusion distance at recombination [49], 


da R‘^ + ^{1 + R) 

0 a^aTUeH 6(1 -|- R^ ’ 


(3.2) 


where Oi-ec is the scale factor at recombination, R is the ratio of energy in baryons to photons, 
rie is the number density of free electrons and ctt is the Thompson cross section. Understanding 
the precise impact of a change in + Aguid is non-trivial [37], since changing H will also 
affect the location of the first acoustic peak, which is extremely well measured. Instead, we must 
simultaneously vary H and Hq, while keeping the angular scale of the first acoustic peak fixed. 
Nevertheless, the physical origin of the effect on the damping tail is still given by (3.2). This 
allows us to understand degeneracies with other cosmological parameters. 

As pointed out in [17], there is an important degeneracy between H and the helium fraction Yp. 
Since helium has a much larger binding energy than hydrogen, increasing (decreasing) the helium 
fraction will decrease (increase) the number of free electrons at the time of recombination, ng oc 
(1 — Yp). It is then possible to change Yp and H simultaneously in a way that keeps the damping 
scale fixed, cf. oc {neH)~^. We therefore expect that the CMB temperature constraints on 

We prefer the use of MCMC techniques over Fisher matrix forecasts because Fisher matrices can underestimate 
the impact of degeneracies on the posterior distributions [47]. We believe this to be the origin of the (small) 
differences between our results and those of ref. [48]. 
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A^eff + -^fiuid and Yp to weaken considerably if we allow both of these parameters to vary. Of 
course, we also must break the degeneracy between A^eff and A^fluid) which are not distinguished 
by their effects on the damping tail. 

Fortunately, future data sets will be much less sensitive to these degeneracies for two reasons: 

• First, as we show in Appendix B, the amplitude of the polarization of the CMB, Qp/, 
is proportional to n~^, but not H, and, therefore, it is sensitive to Yp alone. The key 
feature is that polarization is a direct measurement of the quadrupole at the surface of 
last-scattering, while the damping tail of the temperature spectrum is the integrated effect 
of the quadrupole on the monopole. This difference allows us to break the degeneracy 
between Yp and A^eff -|- Aifluid. 

• Second, as we demonstrated in Section 2, the CMB is sensitive to the perturbations in 
the free-streaming particles and not just their contribution to the background evolution. 
This is illustrated in Fig. 6 and will be explored in more detail in the next subsection. 
What is important here is that the phase shift associated with free-streaming particles is 
not expected to be degenerate with other effects and it is measurable out to very high 
multipoles. Furthermore, as we discuss in Appendix B, the same phase shift appears in 
both temperature and polarization which means there is room for significant improvement 
in the sensitivity to this effect. This will largely eliminate the issues of degeneracies for Agfr- 

3.2 Detecting Free-Streaming Particles 

As we move towards more sensitive experiments, the perturbations in the radiation density will 
play an increasingly important role. In this section, we will demonstrate the potential of current 
and future experiments to detect the free-streaming nature of relativistic species. 

3.2.1 Methodology 

The following data were used in our analysis: 

• The Planck likelihoods were described in detail in [43]. The low-t" likelihoods {2 < i < 29) 
include both the temperature and polarization data (even in the cases labeled “TT-only”). 
For the remaining multipoles {i > 30), we use the Plik joint TT-|-TE-|-EE likelihood. This 
contains information about the TT spectrum up to t'max = 2508 and about the TE and 
EE spectra up to ^max = 1996. The leasing potential is reconstructed in the multipole 
range 40 < i < 400 using SMICA temperature and polarization maps. Eor the TT-only 
constraints, we use the Plik TT-only likelihood with range ‘iO < i < 2508 and the leasing 
reconstruction involves only the temperature map. 

• Our forecasts for CMB-S4 experiments assumed 10® polarization-sensitive detectors with 
a 1 arcminute beam and sky coverage /sky = 0.75. The default observing time was chosen 
to be five years, resulting in a sensitivity of ap = cfpI^/2 = 0.558 /iKarcmin matching the 
most optimistic experimental setup studied in [48]. Our analysis included multipoles up 
to I'max = 5000 for both temperature and polarization. We also studied how our results 
change if we vary the beam size and the maximum multipole. 


19 




V 


EE 

i 


AVf^ 


1950 2000 2050 2100 2150 2200 1800 1850 1900 

i 


1950 2000 2050 


Figure 6. To'p: TT spectrum (left) and EE spectrum {right) for {Nes = 3.046, = 0) 

(blue) and (fVeff = 2.046, A^suid = 1-0) (red) with Vi = i{i + l)C'f/(27r) in units of /rK^. The TT and 
EE spectra represented by the red curves were rescaled by the same constant factor chosen such that the 
height of the seventh peak of the TT spectrum matches for the red and blue curves. The solid and dashed 
lines show the unlensed and lensed data, respectively. The phase shift from and the peak smearing 
from lensing can be seen in both the TT and EE spectra. Bottom: Illustration of the relative difference 
AVg = 8'Diil'Di between the (A^eS = 3.046, A^^auid = 0) and (A^eff = 2.046, A^fluid = 1-0) spectra of the upper 
panels. The green solid and dashed lines are the differences in the unlensed and lensed data, respectively. 
We see that the change is largest in the unlensed EE spectrum. 


Our modification of the Boltzmann code CLASS [50] includes an additional relativistic fluid, 
whose energy density is measured by the parameter Wfiuid defined in (1.2). The equation of 
state and the sound speed of the fluid were fixed to wy = Cy = ^, with initial conditions that 
were chosen to be adiabatic. With this choice, our analytic results imply that this fluid does 
not contribute to the phase shift in the acoustic peaks. We use Monte Python [51] to derive 
constraints on the parameters A^eff and Wfiuid. Whenever the primordial helium abundance Yp 
was not varied independently (which we will refer to as “1^ fixed”), it was set to be consistent 
with the predictions of BBN, using the total relativistic energy density including both A^efi and 
A^fluid in determining the expansion rate. All chains were run until the variation in their means 
was small relative to the standard deviation (using R — l< 0.01 in the Gelman-Rubin criterion). 

Our analysis makes use of the effects of gravitational lensing of the CMB in two distinct ways. 
“Lensing reconstruction” will refer to a reconstruction of the power spectrum of the lensing 
potential from the measurements of the temperature and polarization four-point functions. In 
the case of the CMB-S4 forecasts, the power spectrum of the lensing potential was computed with 
CLASS. CMB lensing also modifies the observed CMB power spectra (TT, TE, EE), primarily 
in the form of smearing the peaks [52], as is illustrated in Eig. 6. “Delensing” removes the effect 
of lensing on these power spectra using the reconstructed lensing potential. This is trivially 
implemented in forecasts in the limit of perfect delensing (we will just output spectra without 
computing the lensing), but is an involved procedure to implement on real data. The utility of this 
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procedure is that lensing moves information from the power spectra to higher-point functions, but 
delensing moves this information back to the power spectra, so that it can easily be accounted 
for in our likelihood analysis (rather than through some more elaborate multi-point function 
likelihood). For a detailed review of gravitational lensing of the CMB we refer to [53]. 

Our Planck 2015 results use the publicly available lensing reconstruction likelihood, but do 
not include any delensing of the power spectra. For CMB-S4, the lensing reconstruction noise 
was computed using the iterated delensing method described in [54] (based on [55-57]). Forecasts 
using delensed spectra assumed perfect delensing, which is a good approximation for a CMB-S4 
experiment across a wide range of multipoles. Taking lensing reconstruction noise into account 
when computing delensed spectra would require a more careful analysis. 

The fiducial cosmology used for all forecasts is described by the following parameters: = 

0.022, = 0.120, h = 0.67, = 2.42 x 10-^, = 0.965, r = 0.078, and A^eff = 3.046. 

3.2.2 Planck 2015 Results 

The Planck 2015 results have reached an important threshold. The level of sensitivity is now 
sufficient to detect the free-streaming nature of the neutrinos (or any additional dark radiation). 
In Table 1, we present marginalized constraints on A^eff and Aiauid; with their posterior distribu¬ 
tions shown in Fig. 7. We show results both for the combined TT, TE, and EE likelihoods, and 
for TT alone. The two-dimensional joint constraints are presented in Fig. 8. In each case, we 
compare the results for fixed Yp with those when Yp is allowed to vary. These results robustly 
demonstrate that there is very little degeneracy between A^efi, A^fluid) and Yp when using both 
temperature and polarization data from Planck. 



TT, TE, 

EE 

TT-only 


varying Yp 

fixed Yp 

varying Yp 

fixed Yp 

A^eff 

^fluid 

2-68l°;i 

< 0.64 

< 0.67 

2 80+0.49 

< 1.08 

9 07 + 0.45 
-0.37 

< 0.94 


Table 1. Best-fit values and la errors for A^eff and 2cr upper limits for ATiuid for the Planck 2015 data. 
Both A^fluid and A"eff are allowed to vary in all cases. The lensing reconstruction and low-P likelihoods 
were used for all of the constraints. 


From the left panels in Figures 7 and 8 we see that the constraints on Ngg are largely insensitive 
to the marginalization over Yp and/or A/juid, even when the polarization data is removed. Since 
A^eff is degenerate with Yp and Aifjuid in the damping tail, this robustness of the constraints 
suggests that they are driven by the phase shift. Since is the unique parameter capable of 
producing the phase shift, the measurement of the latter breaks the degeneracy between Ngg and 
both A^auid and Yp. We also see that adding polarization data leads to a large improvement in 
the constraint on A^eff, most likely because the peaks of the E-mode spectrum are sharper, which 
makes the phase shift easier to measure [17]. This is illustrated in Fig. 6, where we show the 
relative differences in the TT and EE spectra when varying A^eflf and A^uid- While the phase shift 
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Figure 7. Left: Posterior distributions for A^eff from Planck TT, TE, and EE marginalized over A^fluid- 
The black curve involves the marginalization over Yp, while the blue curve keeps Yp fixed. Both likelihoods 
rule out A^eff = 0 at high significance. Right: Posterior distributions for A^fluid from Planck TT, TE, and 
EE marginalized over A^eff- The black curve involves the marginalization over Yp, while the blue curve 
keeps Yp fixed. In both panels, the likelihoods for Planck TT-only with the same marginalizations are 
shown as dashed lines. 
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Figure 8. Left: Constraints on A"eff and A^suid using the Planck TT, TE, and EE likelihoods for varying 
Yp (red) and when Yp is fixed (green). Shown are also the Planck TT-only results for varying Yp (indigo) 
and when Yp is fixed (olive). Right: Constraints on A^eff and Yp using the Planck TT, TE, and EE 
likelihoods for varying A^fluid (red) and when A^fluid is fixed (cyan). Shown are also the Planck TT-only 
results for varying A^fluid (indigo) and when A^fluid is fixed (gray). In both panels, the lensing reconstruction 
and low-P likelihoods were used for all of the constraints. 


is visible in both cases, the size of the effect is larger in the polarization spectrum which increases 
the impact of the E-mode data. 

Similar analyses were performed in [35, 36] using WMAP data (and external datasets). Their 
results are qualitatively similar to our TT-only analysis with Yp fixed, although with weaker 
constraints on N^s and Afluid- By comparison, adding E-mode data further improves constraints 
in the Aeff-A^fluid plane, also when Yp is allowed to vary. 


22 










Figure 8 shows that there is a significant difference in the constraints in the plane, 

with and without the polarization data. Without polarization, allowing iVfluid to vary weakens 
constraints on Ngg and Yp. This is consistent with our discussion of degeneracies in §3.1. Using 
only the temperature data, Yp and Aifjuid are measured mostly from the damping tail, but their 
effects on the latter are degenerate. As a result, the constraints on Yp and weaken when 
A^fluid is allowed to vary. The situation changes when polarization data is added. Now, there is 
very little difference in the constraints as we vary the marginalization over additional parameters. 
This is most noticeable in the right panel in Fig. 8, where the constraints on Yp and Aeff become 
nearly independent of the treatment of Afluid- This feature was anticipated in §3.1, where it was 
observed that polarization breaks the degeneracy between Yp and Aguid. 

In [38], a constraint was recently placed on the effective number of free-streaming species by 
isolating the phase shift in the Planck 2013 temperature data: Agff = 2.3 ^q 4 (68%C.L.) while 
keeping the damping tail fixed and Aeff = 3.5 ± 0.65 (68%C.L.) when marginalizing over the 
effect on the damping tail [58]. To compare to that analysis, we remove the polarization data. 
We then find N^g = 2.89lQg2) which is quite similar to the direct measurement of the phase 
shift.When we add the TE and EE likelihoods, our constraint improves by about a factor 
of two to Aeff = 2.68i*ig33. Erom the estimate^^ in (2.68), we conclude that these constraints 
correspond to a phase shift of about 6 £ ^ 1. While this is compatible with expectations,^^ it is 
nonetheless impressive that the data is sensitive to these small and subtle effects. 

We note in passing that our results show that Yp is somewhat higher than its expected value. 
This is evident in Eig. 8, where the value of Yp is almost Icr higher than the BBN prediction 
of 1), ~ 0.247. This was also observed by the Planck collaboration [5], although their values of 
Yp are somewhat lower.It remains unclear what precisely in the data is driving Yp to larger 
values. This should serve as a reminder that as we increase precision, we will only be able to 
trust a BSM interpretation of our data if accuracy is addressed at the same level. 

3.2.3 CMB Stage IV Forecasts 

While current data is already sensitive to the free-streaming nature of neutrinos, future experi¬ 
ments are expected to improve these constraints by at least an order of magnitude. As we have 
emphasized in the introduction, increasing the sensitivity to the cj(Aeff) ~ 10“^ level probes a 

^^Due to the marginalization over Nnuid and Yp, our analysis is most comparable to the marginalized result: 
Aeff = 3.5 ± 0.65 [58]. One difference in our approach is that it includes information in the amplitude shift 
produced by the free-streaming species. However, this effect is likely subdominant to the phase shift due to the 
degeneracy with the amplitude of the primordial power spectrum. As we will discuss in the next subsection, when 
we allow Yp to vary, future CMB missions get the majority of the sensitivity of Aeff from the phase shift and so 
we expect our methods to produce increasingly similar results. 

^^We can also estimate this more schematically from the knowledge that Aeff = 3.046 produces S£ « 10 for 
£ < 3000 relative to Aeff = 1 [38]. Current constraints allow for roughly a 10 percent variation in Aeff, which 
would imply Si ~ 1. Direct measurements of the individual peak locations are given in [43] with a similar level of 
precision. 

^^Forecasts using the isolated phase shift alone give (j(Aeff) = 0.41 for Planck with polarization [38]. 

^®We found results closer to the Planck values when using CosmoMC/CAMB [59, 60] rather than Monte 
Python/CLASS in the analysis with varying Aeff and Yp. There are several subtle differences in the implemen¬ 
tations of these codes which could be responsible for this disagreement. 
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Experiment 

Delensing 

Reconstruction 

<T{Nes) 

^yp) 

-^fluid 

Planck 2015 

No 

Yes 

0.31 

0.019 

< 0.64 


No 

No 

0.062 

0.0053 

< 0.18 

CMB-S4 

Yes 

No 

0.054 

0.0044 

< 0.17 


Yes 

Yes 

0.050 

0.0043 

< 0.16 


Table 2. Marginalized likelihoods for Planck and CMB-S4. All likelihoods allow Nes, Na^id, and Yp to 
vary. The Planck likelihood uses the TT, TE, and EE power spectra, as well as the leasing reconstruction 
likelihood. The forecasts for CMB-S4 include the possibilities of delensing and leasing reconstruction. We 
excluded the forecast using leased spectra combined with leasing reconstruction since this is known to 
produce overly optimistic error forecasts due to a double counting of leasing information [61]. The double 
counting can be safely ignored for Planck, but will become more important for future experiments [62]. 
Displayed are la error bars for A"eff and Yp, and 2a upper limits for A^fluid- 
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Figure 9. Left: Eorecasted CMB-S4 constraints on Ngg and A^suid for Yp hxed (orange) and varying (blue). 
Right: Eorecasted CMB-S4 constraints on Nes and Yp for A^fluid = 0 hxed (purple) and varying (blue). 



number of plausible BSM scenarios that are currently unconstrained. Preliminary forecasts [48] 
suggest that this level is indeed achievable. Our goal is to extend these results in two ways: (i) 
to include Yp and Afiuid to identify more clearly the types of BSM physics we might be sensitive 
to and (ii) to perform a full likelihood analysis (rather than a Fisher forecast) to ensure that 
degeneracies are treated correctly (see [47] for a discussion). 

The results of our forecasts are summarized in Table 2. Given the constraints from Planck, it is 
not surprising that A^eff is easily distinguished from Afluid with CMB-S4 experiments. As before, 
the constraints on A^eff are significantly stronger than those on Afiuid; which is consistent with the 
interpretation that these parameters are being distinguished by differences in the perturbations 
for the two types of radiation. When both radiation components are included, the detailed 
matching of the acoustic peaks is very sensitive to the phase shifts due to Weflf) but is much 
less affected by A^fiuid- This is illustrated by the left panel in Fig. 9, which presents the joint 
constraints on Wfiuid and A^eff) both for fixed Yp and when it is allowed to vary. When Yp is 
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fixed, there is a strong degeneracy between A^eff and A^fluid which is absent when Yp is allowed to 
vary. Yet, in both cases, the contours close at roughly the same value of A"eff, indicating that the 
degeneracy is broken in a way that is insensitive to the damping tail. 

As we have emphasized throughout, the strong constraint on A^eff arises in part from the 
sharpness of the peaks of the E-mode spectrum, which leads to better measurements of the phase 
shift. This same intuition explains why the constraints on Ngg are strengthened by using the 
delensed power spectrum. One of the well-known effects of lensing is a smearing of the acoustic 
peaks [52] (see Fig. 6). By delensing the power spectra, we sharpen the peaks, which makes the 
measurements of the phase shift more precise. We illustrate the impact that delensing can have in 
Fig. 6, where we see that the effect of changing Agfr produces a much larger relative change on the 
unlensed data. This procedure is analogous to reconstructing the BAO peak to sharpen distance 
measurements [63]. In both cases, non-linearities transfer information from the power spectrum 
to higher-point correlation functions. Delensing or BAO reconstruction moves that information 
back to the power spectrum, so that it is more easily accounted for in these analyses. As a result, 
the error in the measurement of any quantity that is sensitive to the sharpness of the peak will 
be reduced (like the phase shift and the BAO scale). This suggests that delensing will be a useful 
tool for improving constraints on cosmological parameters in these future experiments. 

While much of our efforts have been devoted to understanding the degeneracies between Ngg, 
A^fluid, and Yp, the actual physical models we wish to constrain may not exhibit these degeneracies. 
First of all, many models with additional light fields still have A^fluid = 0 (i.e. only free-streaming 
radiation) [25], which would in principle allow us to combine information from the damping tail 
and the phase shift to constrain A^eff- Furthermore, while Yp is often affected by BSM physics 
at the time of BBN, the precise degeneracy needed to keep the damping tail fixed is unlikely to 
occur naturally. The constraints on such models would be significantly stronger. As we see in 
Table 3, a factor of 3 to 4 improvement in the constraints is possible when this degeneracy with 
changes in Yp is not present. 


Experiment 

^{Nes) 



0.30 

0.018 

Planck 2015 


0.19 

— 


0.048 

0.0027 

CMB-S4 


0.013 

— 


Table 3. Results for the marginalized Icr errors for Planck (TT, TE, EE, lensing reconstruction) and 
forecasts for CMB-S4 for Nes and Yp, when Yfluid = 0 is held fixed. The CMB-S4 forecasts assumed both 
delensing and lensing reconstruction. A dash in the cr(Yp) entry indicates that Yp was fixed by consistency 
with BBN. 


Finally, perhaps the most important feature of these forecasts is the dramatic improvement, 
relative to Planck 2015, that can be expected from a plausible experimental configuration of CMB- 
S4. Our projections suggest that a factor of 5 to 10 improvements are achievable, but we should 
also investigate the robustness of this conclusion to changes of the experimental configuration. 
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Parameter 

1' 

2' 

3' 

W = 3000 

W = 4000 

o-(Aefr) (Tp fixed, Nff^^ff = 0) 

0.013 

0.015 

0.016 

0.023 

0.015 

o-(Aefr) (Tp fixed, Afiuid / 0) 

0.026 

0.027 

0.029 

0.034 

0.028 

o-(Aefr) (Tp varying, Nff^id = 0) 

0.048 

0.051 

0.055 

0.058 

0.052 

cr(Aefr) (Tp varying, A^fiuid / 0) 

0.050 

0.052 

0.055 

0.061 

0.051 

A^fluid (Tp varying) 

< 0.16 

< 0.17 

< 0.18 

< 0.20 

< 0.17 

A^fluid {Yp fixed) 

< 0.068 

< 0.072 

< 0.076 

< 0.090 

< 0.072 


Table 4. Forecasts for a CMB-S4 experiment with varying beam size and maximum multipole assuming 
10® detectors. We take £max = 5000 as we vary the beam size and 9b = V when we vary .^max- Displayed 
are Icr error bars for A^eff and 2a upper limits for iVfluid- 


Here, we study variations of the beam size {9f,) and the maximal available multipole (-^max)- These 
are important for two reasons: (i) the beam size is nltimately a choice made within the context of 
limited resonrces and (m) the presence of foregrounds or systematics make it difficult to predict 
t? m ax reliably beforehand. In Table 4 we show the forecasts for various values of 6b and .^max, 
assuming 10® detectors. 

There are two simple lessons we can draw from this table. First of all, measuring modes 
with i > 3000 offers only a moderate improvement on onr constraints. Similarly, the benefit to 
redncing the beam size significantly is also limited. Most of the improvement in the measurement 
of the phase shift is coming from high-precision measurements of E-modes with i < 3000. In 
contrast, when Yp is fixed and Wfiuid = 0) the sensitivity to the beam size and I’max is stronger, 
which suggests^® that we are gaining nseful information from the damping tail at £ > 3000. 

We have fixed the number of detectors to be 10®, but the precise value will play a very 
important role in the ultimate reach of CMB-S4, not just for N^ff, but for most physics targets. 
For Neff specifically, we are still far from the limit set by cosmic variance in the E-mode power 
spectrum for i > 1000, which is where most of the improvements in the phase shift measurements 
are coming from. As a result, within the range of 10“^ to 10® detectors being considered, we 
improve constraints significantly by further reducing the detector noise and, hence, increasing 
the number of detectors. Eorecasts for Ngs with varying numbers of detectors were studied 
in [48], which confirm this intuition. 

3.3 Time Evolution of Radiation Densities 

One of the important benefits of measnring the fluctuations in the dark radiation is that it 
eliminates the degeneracy between Yp and Neff. This has important consequences for constraints 
on BSM physics, since Yp is sensitive to the total radiation density at the time of BBN (3 min 
after the Big Bang), while Neff and iVfluid are related to these radiation densities around the time 
of recombination (380 000 years after the Big Bang). As a result, CMB measurements of Yp and 

^®We have not been careful to account for foregrounds in temperature or polarization. Information in our 
forecasts coming from high-£ temperature data is unlikely to be available in a real experiment. 
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A^eff probe scenarios where these densities change between those two times, perhaps through the 
decay of a heavy particle or through some other production mechanism [64]. 

It is useful to translate constraints on Yp into constraints on the radiation density at the time 
of BBN [65], 

« 0.247 + 0.014 X <5iVeffXd , (3-3) 

where — 3.046. Marginalizing over Alguid, we found an error in the 

helium fraction of criYp) = 0.0043 for CMB-S4 (see Table 2), while for fixed A^fiuid = 0, we get 
cr{Yp) = 0.0027 (cf. Table 3). Using (3.3), these constraints imply 




0.31 iVfluid + 0 
0.19 iVfluid = 0 


(3.4) 


The last constraint is stronger than the current best limit of = 0.28 from BBN 

alone [4]. The CMB will therefore provide independent measurements of Neff at two different 
times in a single experiment, each surpassing our current level of sensitivities from combining 
multiple probes. 

While the constraint on from a CMB-S4 experiment is only a modest improvement 

over current measurements from primordial abundances, it has the unique advantage that it is 
a clean measurement (i.e. it is not affected by astrophysical processes at later times) and it can 
be combined with measurements of other cosmological parameters (e.g. VL},h?) without combining 
different data sets. A common approach with current data is to combine the constraints from 
the CMB and primordial abundances in order to improve the overall sensitivity to Neg, in the 
case where it is time independent. From Table 3, we see that if we do not allow a variation 
in Neff between BBN and the CMB, we get very strong constraints on Neff due to the lack of 
degeneracies in the damping tail. These results are sufficiently strong so that it is unlikely that 
including information from primordial abundances will lead to much improvement. 

When discussing the time variation of Neff, we have focused only on the model-independent 
measurement implied by varying Neff and Yp independently. As a result, the constraints we derive 
are controlled primarily by the degeneracy between Neff and Yp in the damping tail. Without 
this degeneracy, the constraints on Ngff are much stronger. In realistic models, it may be the case 
that both Neff and Yp are changed independently, but that they do not produce this degeneracy 
in the relevant range of parameters. For such models, a dedicated analysis of CMB data would 
likely offer a much larger gain over the current limits from primordial abundances. 


3.4 Implications for BSM Physics 

In this section, we have focused on simple descriptions of BSM physics in terms of the effective 
parameters Neff, Nffuid, and Yp. These parameters capture important aspects of CMB physics 
and our cosmological history, which may be used to test a number of scenarios for BSM physics 
(see Table 5): 


• For minimal extensions of the Standard Model with a light field, e.g. [25, 30, 42], Ngg is 
sensitive to the freeze-out of the particle. At current levels of sensitivity, a{N^s) ^ 0.1, we 
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Signature Influenced by Degeneracies broken by 


CMB damping tail 


Phase shift, Polarization 

Phase shift 

iVeff, iSFf 

Scale dependence 

Spectral distortions 

ppost-BBN 

-^inj 


Primordial abundances 

afBBN _i_ afBBN ]\t ^BBN tt'BBN 

CMB 


Table 5. Cosmological probes of BSM physics and their sensitivity to free-streaming and non-free- 
streaming radiation {NeS and A^suid), the number of active neutrinos (iV^), the baryon-to-photon ratio ( 77 ), 
and the amount of energy injection (-Einj)- The superscripts BBN, CMB, or post-BBN denote the time at 
which a quantity is being probed, where post-BBN refers to redshifts of z < 10® when spectral distortions 
become possible. The parameter abstractly stands for isocurvature fluctuations, while jSTf denotes 

modified gravity. 

can rule out some scenarios where particles freeze-out after the QCD phase transition [25]. 
Freeze-out before the QCD phase transition typically dilutes the contribution to Ngg by 
a factor of 10, which allows such models to easily evade current constraints. Fortunately, 
these scenarios are likely to be accessible with CMB-S4 experiments [30]. For these cases, 
we are sensitive to sufficiently early times so that BSM physics above the TeV scale may 
be important. 

• Measurements of the effective number of free-streaming particles at recombination, 

are also sensitive to energy injected into the Standard Model particles after the time of 
neutrino decoupling (~ 0.1 sec). Depending on the time and nature of this energy injection, 
it may alter the primordial abundances or introduce spectral distortions which would dis¬ 
tinguish it from a new light field. For example, a decay to photons after BBN would lower 
and 77 ®^^ (the baryon-to-photon ratio at BBN), while keeping the radiation density 
at BBN, fixed [22]. 

• Energy injection of many kinds is a typical byproduct of changing iVeff, but may also be the 
dominant signature of BSM physics. Decays during BBN can disrupt the formation of nuclei 
without substantially changing the total energy in radiation. Alternatively, recombination 
is very sensitive to energy injection [ 6 ] which can alter the form of the visibility function. 

• As we discussed in Section 2 , phase shifts of the acoustic peaks may also be produced by 

isocurvature perturbations (denoted by in the table). We offered a simple curvaton- 

like example of this effect, but we expect to be broadly sensitive to physics in the dark sector 
that is not purely adiabatic. Since there are many good reasons to imagine why isocurvature 
perturbations might arise in the dark sector, this motivates a future exploration of the 
observability of these effects. 

• Finally, we have assumed the validity of the Einstein equations throughout. This enforced 
that d>_ = 0 in the absence of anisotropic stress. Modified theories of gravity are often 

^^The common element of both of these examples is that the tail of the Boltzmann distribution is playing a 
critical role (due to the large value of ?7~^). As a result, the change to the small number of high-energy photons is 
more important than the total energy density. 
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parameterized in terms of their change to the Einstein constraint equation and the corre¬ 
sponding effect on <h_; see e.g. [66]. Since the field <h_ played an important role in our 
analysis of the phase shift, it would be interesting to explore how the result changes for 
specific modifications of GR. Conversely, the phase shift of the CMB spectrum may be an 
interesting probe of modified gravity. 

4 Conclusions and Outlook 

CMB observations have become precise enough to probe the gravitational imprints of BSM physics 
on the perturbations of the primordial plasma. In the upcoming era of CMB polarization ex¬ 
periments, our sensitivity to these subtle effects will increase significantly and will offer new 
opportunities in the search for new physics. It is therefore timely to re-evaluate how CMB data 
can inform our view of the laws of physics. 

In this paper, we have explored how the phase shift of the acoustic peaks might be used as such 
a probe. This phase shift is particularly interesting because analytic properties of the Green’s 
function of the gravitational potential strongly limit the possible origins of such a shift to 

i. waves propagating faster than the sound speed of the photon-baryon fluid, 
ii. isocurvature fluctuations. 

For adiabatic initial conditions, the phase shift is most easily generated by free-streaming radi¬ 
ation and becomes an excellent measure of the effective number of neutrinos iVeff at the time of 
recombination. Realistic models of isocurvature fluctuations typically produce a scale-dependent 
phase shift, which allows them to be distinguished from changes to the energy density of the 
radiation. 

What makes these results particularly compelling is that current and future CMB experiments 
are sensitive enough to detect these phase shifts at high significance [38]. We have demonstrated 
this with an analysis of the 2015 data from the Planck satellite and forecasts for a CMB Stage IV 
experiment (see Fig. 10). Our results include a clear detection of the free-streaming nature of 
neutrinos. They also highlight the important role played by the polarization data in breaking the 
degeneracy between the contributions from free-streaming and non-free-streaming species, Vgff 
and Vfluid, as well as that with the helium fraction Yp. 

The increased sensitivity to differences in the evolution of perturbations motivates revisiting 
the predictions of specific models of BSM physics. Reference [25] performed a comprehensive 
study of ‘minimal’ effective field theory realizations of technically natural forms of dark radiation. 
However, the predictions of these models were only compared to existing constraints on Ngg, 
i.e. only modifications to the background density were taken into account. (Recently, these 
estimates were updated [30] in light of the higher sensitivity expected in CMB-S4 experiments.) 
In the future, we plan to revisit this approach, including the effects of the perturbations. These 
effects may arise from new degrees of freedom (e.g. [22-30]) or from non-standard properties 
of the SM neutrinos (e.g. [39-41]). We also hope to explore ‘non-minimal’ extensions of the 
effective theories of [25]. For example, decays of massive fields are abundant in well-motivated 
extensions of the Standard Model, but were forbidden in [25] on the basis of being non-minimal. 
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Figure 10. Left: Planck constraints on the effective number of free-streaming and non-free-streaming rel¬ 
ativistic species, TVeff and fVfluid, allowing the helium fraction Yp to vary (red) and keeping it fixed (green). 
Right: Planck constraints on Yp and Neff for varying iVfluid (red) and when Nfluid = 0 is fixed (cyan). In 
both plots, the contours from the CMB-S4 forecasts, presented in Fig. 9, have been included to show the 
expected improvements in sensitivity. 

As we illustrated, the type of energy injection that would be produced from such a decay can 
be observable over a wide range of times, even in the absence of new light fields. Fnrthermore, 
the detailed predictions for A^eff may be significantly altered from the minimal case, which can 
change the impact of future constraints. 

The expected improvement in the constraints on extra relativistic species in future experiments 
is quite remarkable, e.g. a{Nes) ~ 0.01. This provides an opportunity to probe BSM physics at 
a much more precise level than was previously possible. We are optimistic that this will teach 
us something interesting. We will either discover a whole new world of dark physics, or learn to 
what remarkable degree it is decoupled from the rest of physics. 
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A Comments on Matter 


In the main text, we computed the phase shift of the photon density fluctuations assuming a 
radiation-dominated background. While this simplihcation made an analytic treatment possible, 
we may wonder if it misses important effects. In this appendix, we will bridge this gap to the 
degree which is possible without using numerics, focusing on the contributions from free-streaming 
radiation. 

There are several reasons why we want to understand the contributions to the phase shift from 
modes in the matter era. First, recombination occurs during matter domination and, therefore, 
in principle, it could be important for every mode in the CMB. Second, modes corresponding 
to large angular scales (small €) enter the horizon during (or near) matter domination and their 
complete evolution is therefore governed by the physics in the matter era. Finally, ref. [38] found 
a logarithmic dependence of the phase shift on i for observable modes. One may be tempted to 
interpret this as an effect of the finite matter density. Our goal in this section is to further clarify 
these effects, by studying the limits I ^ oo and £ —)> 0, accounting for the contributions from 
matter. We will study these limits in turn: 


• We first consider modes which entered the horizon during the radiation era. These cor¬ 
respond to small angular scales in the CMB anisotropy spectra. We begin by writing 
eq. (2.20) as 


rVeq ry 

B{y)= dy $+(?/') cosy -b / dy'^>+(y') cosy' , 
Jo Jyec{ 


— -Srad 


= B 


mat 


(AT) 


where yeq = c^kTeq corresponds to the moment of matter-radiation equality. Modes that 
entered the horizon long before Teq correspond to yeq 1. For these modes, the first term 
in (A.l) can be approximated as 


1*00 

-Brad- / dy'^>Y'"‘^^(y') cosy', (A.2) 

Jo 

which is precisely the result computed in Section 2. The main correction from the matter 
era then is the second term in (A.l): 


Bmat ^ f dy' (y') COS y'. 

yeq 


(A.3) 


To estimate this effect, we simply have to repeat the discussion of §2.4.2 for the matter era. 
The important difference is that ex ^ px/p now is not a constant, but scales as a~^ oc t“^. 
Setting Tin in the matter era, we have ex = ex,inand eq. (2.55) becomes 


8k^ 16 1 

^-{y) = -^ex,in2/inC^x(y) = — ^^ex;myfnDx,2{y) ■ (A.4) 

To determine <I>_(y) to first order in ex, we only need the quadrupole moment Dx, 2 {y) to 
zeroth order. From (2.58), we get 

B%{y) = dx;mj2[c~^{y -y-m)] 

+ [ dy' I [c-^{y - y')] - ^j3 [c-^{y - y')] | , (A.5) 

W '^yin V J 
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where we have used that = const, during the matter era. In the limit, y 3> 1, this 
leads to 

^ 16ex,iny?n - -7 (^+.in + \dx,m) • (A.6) 

Hence, we get $_ oc y~^ —)• 0 in the limit y —)• oo. At late times, $+ is therefore no longer 
sourced by <h_ and will be given by the homogeneous solution (with coefficients that may 
depend on ex)- Importantly, the value of B will be the same as that predicted in Section 2. 
Hence, at high k (and thus £), the phase shift of the acoustic oscillations will be equal to 
the value in a radiation-dominated universe. 

• Next, let us study modes which entered the horizon during the matter era, corresponding 
to large angular scales in the CMB. In this case, it is more challenging to cleanly separate 
the result into a correction to the amplitude of oscillations and a phase shift. Our primary 
goal will be to understand how the result scales with wavenumber k in the limit A: —?■ 0. 
Fortunately, this scaling is the same for the amplitude correction and the phase shift and 
is easy to understand analytically. 

Intuitively, we expect the contributions from dark radiation (including neutrinos) to vanish 
as A: —)• 0. As we lower k, the time of horizon entry increases compared to the time of matter- 
radiation equality and, therefore, the radiation energy density should be diluted relative to 
the matter. Since this radiation only affects observations through its gravitational influence, 
its role in the evolution of the modes should become negligible. 

We can confirm this intuition by returning to (A.4) and noticing that yin = c^kr^^, where 
Tin is a fixed time which is independent of k, e.g. we may choose Tin to be the time of 
matter-radiation equality. We therefore have = ex,inC^A:^Ti^ y(y), and the correction to 
d-y at linear order in ex will take the form 

= ex,in c^k‘^T^n [ diy'f{y,y'). (A.7) 

dpin 

If the integral converges as y —)■ oo, it is clear that dj^ oc A:^ —)■ 0. In fact, if the integral 
diverges as r —)■ oo, the result will be suppressed by additional powers of k, due to the 
scaling of the upper limit of integration (y = c^kr at fixed r). Hence, we conclude that the 
amplitude and phase corrections from neutrinos (or any dark radiation) will vanish at least 
as fast as A;^. 

From these asymptotic scaling arguments, we draw the following conclusions: 

• For A; —)• oo, the phase shift due to free-streaming particles approaches a constant. 

• For A; —)• 0, the phase and amplitude corrections scale at least as A;^. 

In the flat-sky limit, these results translate approximately to £ ~ A;(ro — Tree), where tq is the 
conformal time today. We therefore expect a constant phase shift at high i. Given that matter- 
radiation equality corresponds to relatively low £, we do not expect our asymptotic formula for 
A: —>• 0 to be more than a rough guide. The primary purpose of this discussion was to highlight 
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that a power law is the likely behavior, simply due to the power law decay of the energy density 
of the extra radiation. As a result, the phase shift per I should be some function that interpolates 
between a power law and a constant and is unlikely to follow the ansatz of [38] in detail (although 
the logarithmic dependence appears to work well enough on intermediate scales). 


B Comments on Polarization 


The analytic discussions in the main text were phrased in terms of the temperature anisotropy, 
but as we saw in Section 3, CMB polarization plays a crucial role in present and future data 
analysis. In this appendix, we show that the phase shift of the polarization spectrum is the same 
as that of the temperature spectrum. 


Following [49], we write the Boltzmann equation for the amplitude of polarized anisotropies. 


0p, as 


Qp + ikfj,Qp = —k 


-0p + -(i-P2(i^))n 


(B.l) 


where 11 = 02 + 0p,o + &p ,2 and k = —neapa is the time derivative of the optical depth k 
(to avoid confusion with the conformal time r). The temperature quadrupole is determined by 
the photon anisotropic stress, 02 = Equation (B.l) admits a solution as a line-of-sight 

integral, 

0p(ro) = y ° dr e^kt,{r-ro)-K{r) (^^k{T) (^2 _ 1) u{t)^ . (B.2) 


The integral in (B.2) is proportional to the visibility function —ke~^ and is, therefore, peaked at 
the surface of last-scattering. In the limit of instantaneous recombination, —ke~^ ~ Spip — Tree), 
we get 

0p(ro) ~ - /i2)n(rrec) • (B.3) 

Solving for Ft to leading order in k 3> 1, one finds 11 ~ |02 — —^kk~^ 0i (using the collision 
term in the Boltzmann equation for temperature). Applying the continuity equation, = 
—3kQi, and performing a multipole expansion, one finds 


©p.K'To) - 


5 • 

— d^{k,Trec) k~^{T,ec) 



je{kTo) . 


(B.4) 


Two facts should be noticed about this result: 


• Qp^i (X d-y .—Since the time derivative will not affect the phase shift from dark radiation, 
we see that the locations of the acoustic peaks in the polarization spectrum are affected by 
the fluctuations in the dark radiation in the same way as in the temperature spectrum. 

• Qp^e (X k~^ (X n~^. —This is important because it allows the degeneracy between H and Ug 

(or Yp) in the damping tail (which scales as cf. §3.1) to be broken. 
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